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Abstract 


The application of an optimization method to 
determine the propeller blade twist distribution 
which maximizes propeller efficiency is presented. 
The optimization employs a previously developed 
method which has been improved to include the ef- 
fects of blade drag, camber and thickness. Before 
the optimization portion of the computer code is 
used, comparisons of calculated propeller efficien- 
cies and power coefficients are made with experi- 
mental data for one NACA propeller at Mach numbers 
in the range of 0.24 to 0.50 and another NACA pro- 
peller at a Mach number of 0.71 to validate the 
propeller aerodynamic analysis portion of the com- 
puter code. Then comparisons of calculated propel- 
ler efficiencies for the optimized and the original 
propellers show the benefits of the optimization 
method in improving propeller performance. This 
method can be applied to the aerodynamic design of 
propellers having straight, swept, or nonplanar 
propeller blades. 


Introduction 


Optimization techniques are receiving in- 
creased attention as a method of solving real world 
problems and are being regarded as a unique tool 
for analysis, design, and selection processes. 

There are wide spread applications of optimization 
methods in science, engineering, and a variety of 
other fields. In this paper, an optimization 
method is applied to maximize the propeller effi- 
ciency at a pre-specif ied design condition. Previ- 
ous investigations for optimizing propel ler twist 
used the method of calculus of variations®*^ and a 
classical method.® All of these approaches are 
formulated for straight propeller blades. The op- 
timization method presented herein is an improved 
version of the method of Ref. 9, where it was dem- 
onstrated by analyzing propellers with tip devices 
that the method could be applied to complicated 
blade shapes. Features not available in the method 
of Ref. 9 but included in the optimization program 
used in this paper are the effects of blade drag, 
camber, and thickness. These additions improve the 
analysis and account for blade profile drag, there- 
by allowing the calculated efficiencies to be com- 
pared with experimental data. 

The objective of this paper is to compare some 
calculated results of the propeller performance 
with the experimental data and then determine the 
maximum propeller efficiency by optimizing the 
blade twist distribution. 


Description of the Propeller Analysis and the 
— ~ Optimization Method 


The optimization computer code includes pro- 
peller performance analysis coupled with optimiza- 
tion procedures. In this paper, the computer code 
was applied to straight blade propellers with uni- 
form inflow velocity field. However, the method 
can handle propellers having an arbitrary blade 
shape (swept or nonplanar propeller blades), any 
number of blades, and operating in a uniform or 
nonuniform flow field. The method currently works 
for a single rotating propeller but is being ex- 
panded to handle counter-rotating propellers. The 
performance analysis is based on the curved lifting 
line method. 9 Summaries of the curved lifting line 
method and its use in the analysis of advanced pro- 
pellers can be found in Refs. 10 to 12. Therefore, 
only a brief description of the curved lifting line 
method and the optimization procedures will be 
given herein. The propeller blade is replaced by a 
bound vortex and control points, which are located 
at the quarter chord point and three quarter chord 
point respectively. The propeller is modeled by 
blade elements, each of which has a blade angle, a 
chord length, a thickness ratio, and a design lift 
coefficient or camber. The varying vortex strength 
along the radius causes vorticity to be shed as 
helical vortex filaments. Influence coefficients 
are calculated, and the circulations are solved by 
satisfying the boundary condition at the control 
points. The lift coefficient of the blade at any 
radial location is then determined from the bound 
vortex strength at that radius. Drag is provided 
by correlated two dimensional airfoil data. 

The optimization technique used was the 1970 
version of SUMT (Sequential Unconstrained 
Minimization Technique) based on Fiacco and 
McCormick's penalty function method. ”,14 A 
variety of optimization methods are available but 
the selection of the penalty function method”, 14 
was due to its high reliability. The penalty func- 
tion of Fiacco and McCormick can be stated as 

i m ? k 

P(x,r) = f(x) + - E h.(x) - r £ In (g^x)), 
r 1-1 1 j-1 J 

where 

P(x,r) penalty function 

x design vector (angle of twist of the 

blade elements) 

f ( x ) objective function (the negative value 

of the propeller efficiency) 
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r penalty parameter, a positive monotoni- 

cally decreasing sequence, r l > r t+1 
h^fx) equality constraints 

gj(x) inequality constraints 

m total number of equality constraints 

k total number of inequality constraints 

The problem is to minimize P(x,r) by solving 
a sequence of unconstrained minimization problems. 
The search direction was calculated by the method 
of Davidon-Fletcher-Powel 1 . 14-16 optimal step 
length is calculated by the golden section search 1 '* 
which locates the current minimum by regional elim- 
inations. A complete description of the optimiza- 
tion method to obtain the best propeller blade 
twist is given in Ref. 9. 

Prope 1 1 er Geometry 

Two NACA propellers which had experimental 
wind tunnel test data available were selected to 
compare the analytical and optimized results. The 
geometries of these two propellers are described 
below. The propeller blade geometries for the two 
propellers are shown in Fig. 1. Both of them were 
built with NACA 16 series airfoil sections. The 
lower speed propeller, ' (Fig. 2) designated by 
NACA 10— ( 3) (08) -03, is a two blade propeller, 
tested at speeds up to Mach number 0.64. Experi- 
mental results shown in Ref. 17 present the per- 
formance of this propeller in terms of apparent 
efficiency. However, because of the relatively 
low speed of the test data and the small size of 
the body aft the propeller, the difference between 
net and apparent efficiencies is sufficiently small 
to warrant the comparison to calculated net effi- 
ciency made herein. 

The higher speed propeller, 18 (Fig. 3) desig- 
nated by NACA 4— ( 3 ) ( 08 ) —045 , was tested for speeds 
up to Mach 0.73. It is different from the NACA 
10-( 3 ) (08) -03 propeller mainly in the design lift 
coefficient and the chord length distribution. 

For reference purposes, the NACA propeller 
designation numbering system (e.g., NACA 
4-(3)(08)-045) is explained below. The number of 
the first group is the diameter in feet; the number 
of the second group (enclosed within the first set 
of parentheses) is the design lift coefficient (in 
tenths) of the blade section at the 0.7-radius 
station; the number of the third group (enclosed 
within the second set of parentheses) is the thick- 
ness ratio of the blade section at the 0.7-radius 
station; and the number of the fourth group is the 
blade solidity expressed as the ratio of the blade 
chord at the 0.7-radius station to the circumfer- 
ence of the circle having a radius 0.7 of the 
propeller tip radius. The NACA 4- (3) (08) -045 pro- 
peller thus has a diameter of 4 ft and the blade 
section at the 0.7-radius station has a design lift 
coefficient of 0.3, a thickness ratio of 0.08, and 
a blade solidity of 0.045. 

Aerodynamic Performance 

The computer code was first used, without the 
optimization procedure, to calculate the aero- 
dynamic performance of the two propellers using 
the original twist distributions. The geometry of 
the lifting line used in the analysis was calcu- 
lated using the geometry of the original design as 


shown in Fig. 2 for the lower speed propeller and 
Fig. 3 for the higher speed propeller. 

For the lower speed propeller, the blade angle 
at 0.75 radius was specified in Ref. 17 tobe 35* 
for Mach numbers from 0.24 to 0.35, and 45° for 
Mach numbers from 0.35 to 0.50. These Mach numbers 
correspond to advance ratio from 1.2 to 2.6. The 
predicted results are compared with test data 1 ' in 
Figs. 4 and 5. In Fig. 4, the predicted power co- 
efficient is compared with experimental data as a 
function of advance ratio. The comparison shows 
that both the shape and the level of the power co- 
efficient are predicted very well. Differences are 
within one percent. Figure 5 presents the pre- 
dicted and measured efficiencies. Again, the com- 
parison shows that both the shape of the efficiency 
curve and the efficiency level agree well. The 
maximum difference is 1.5 percent. 

The performance of the higher speed propeller 
was calculated at different values of advance ratio 
from 3.1 to 3.8 for a blade angle at 0.75 blade 
radius of 55°. The test Mach number for the higher 
speed propeller 18 was not corrected for tunnel wall 
interference. When this correction is made, the 
tunnel datum Mach number of 0.70 corresponds to a 
freestream Mach number of 0.71. The analytical 
calculations were done at this freestream Mach 
number of 0.71. In Fig. 6, the predicted power 
coefficient and the experimental data are compared. 
For lower values of advance ratio, the predictions 
agreed well with the data, but for advance ratios 
above 3.4 the predictions show a slight overesti- 
mate for power coefficient. A comparison of the 
predicted and experimental efficiency values for 
the higher speed propeller is shown in Fig. 7. 

The predicted efficiency underestimates the data 
by about one percent for advance ratios below 3.5, 
and overpredicts the efficiency at advance ratios 
above 3.5. 

All these comparisons between the analytically 
predicted results and the test data demonstrate the 
capability of the code for performance analysis and 
validate its use for the optimization. 

Results of the Optimization 

Optimization performed on the lower speed and 
the higher speed propellers were based on a typical 
cruise condition for each propeller. For the lower 
speed propeller, the cruise condition was at a 
freestream Mach number of 0.44, an advance ratio 
of 2.2 and a blade angle at the 0.75 radius of 45°. 
For the higher speed propeller, a typical cruise 
condition was at a freestream Mach number of 0.71, 
an advance ratio of 3.3 and a blade angle at the 
0.75 radius of 55°. For each propeller, the ini- 
tial input of the blade twist for the optimization 
process was the twist of the original design. The 
equality constraint for the optimization process 
was the power coefficient at cruise. 

In Fig. 8, the original and optimized propel- 
ler blade twist distribution for the lower speed 
propeller is shown. The shape of the optimized 
twist and original design are almost the same. 

Figure 9 shows the comparison of the propeller 
efficiencies that were calculated for the original 
and optimized lower speed propeller design. To 
show the magnitude of the drag loss, the ideal ef- 
ficiency was calculated by an option in the code 
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and shown on this figure. The ideal efficiency 
accounts for the swirl loss, tip loss, and axial 
momentum loss but does not include the drag loss. 
The blade angle at the 0.75 radius was varied at a 
constant advance ratio so that the value of effi- 
ciency at different values of power coefficient 
could be obtained. As shown in Fig. 9, there is 
actually no improvement in efficiency by the opti- 
mization process. This lack of improvement is not 
surprising since there was almost no change in 
twist. It thus appears that the original propeller 
blade twist design was very good. 

Figure 10 shows the twist distribution for the 
original design and the optimized blade of the 
higher speed propeller. The optimized blade shows 
a different distribution than the original design, 
having higher relative blade angles very near the 
hub and in the outboard regions. Blade angles are 
smaller in the middle part of the optimized blade 
than in the original design. Thus the optimized 
propeller blade would have higher loading at the 
hub and over the outer 20 percent of the blade, 
but a lower loading across the mid-portion of the 
blade. 

A comparison between the predicted efficien- 
cies of the optimized and original blade designs 
is presented in Fig. 11. The comparison was made 
for different values of the power coefficients as 
was done in Fig. 9. The optimization process im- 
proved the higher speed propeller efficiency by 
about 1 percent at the cruise power coefficient of 

0. 18. At lower power coefficients, improvements 
as large as 2 percent were predicted. Again, the 
ideal efficiency was added on the figure for a 
comparison of the blade drag. The larger drag 
losses on the higher speed propeller as compared 
with the lower speed propeller (Fig. 9) are due to 
the higher compressible drag caused by the much 
higher helical tip Mach numbers. 

Conclusions 

A method of maximizing propeller performance 
has been presented and demonstrated using the blade 
geometry of two NACA propellers. The method in- 
cludes the effects of blade drag, camber, and 
thickness. Comparisons made in this paper show 
that the present method can be employed to perform 
propeller aerodynamic design. Also, the optimiza- 
tion method provides a tool for decreasing the 
time it takes to design a propeller. The paper 
presents the results of optimization of straight 
blade propellers but the method is equally capable 
of optimizing swept or other nonplanar propeller 
blades. 
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BLADE WIDTH RATIO, 
b/d, AND BLADE 
THICKNESS RATIO, t/b 


LOWER SPEED (LS) PROPELLER 
(NACA 10-(3) (08) -03) 

(NACA TN 2859) 


HIGHER SPEED (HS) PROPELLER 
(NACA 4-(3) (08)-045) 

(NACA REPORT 999) 


Figure 1. - Propeller blade planforms used in the analysis. 
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Figure 2. - Blade geometry for the lower speed (LS) propeller (NACA TN 2859). 
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Figure 3. - Blade geometry for the higher speed (HS) propeller (NACA report 999). 
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Figure 4 - Comparison of analytical and experimental power 
coefficients for the lower speed propeller. Two blade propeller; 
rotational speed, 1350 rpm. 




Figure 5. - Comparison of analytical and experimental apparent 
efficiencies for the lower speed propeller. Two blade propeller*, 
rotational speed, 1350 rpm. 



Figure 6. - Comparison of analytical and experi- 
mental power coefficients for the higher speed 
propeller. Freestream Mach number, 0.71; 
two blade propeller; blade angle at 0. 75 radius, 
55°. 
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Figure 7. - Comparison of analytical and experimental 
net efficiencies for the higher speed propeller. Free- 
stream Mach number, Q. 71; two blade propeller; 
blade angle at Q. 75 radius, 55°. 



Figure & - Comparison of the original and optimized 
blade-twist for the lower speed propeller. Freestream 
Mach number, 0.44; two blade propeller; advance 
ratio, 2. 2; blade angle at 0. 75 radius, 45°. 
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Figure 9. - Comparison of predicted efficiencies 
for the lower speed propeller. Freestream Mach 
number. 0.44; two blade propeller; advance 
ratio, 2.2; blade angle at 0.75 radius, 45°. 



Figure 10. - Comparison of the original and optimized 
blade-twist for the higher speed propeller. Free- 
stream Mach number, 0.71; two blade propeller; 
advance ratio, 3. 3 ; blade angle at 0. 75 radius. 55°. 
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Figure 11. - Comparison of predicted efficiencies 
for the higher speed propeller. Freestream 
Mach number, 0.71; two blade propeller; ad- 
vance ratio, 3.3; blade angle at 0.75 radius, 55°. 
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